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Abstract 

In this article we study quasilinear systems of two types, in a domain SI of M. N : with 
absorption terms, or mixed terms: 

^ { AqV =S\ W{X = =X 

where S, /i > and 1 < p, q < N, and D = 5/x — (p — l)(q — 1) > 0; the model case is 
A p = A p , A q = A g . Despite of the lack of comparison principle, we prove a priori estimates 
of Keller-Osserman type: 

t>(q — 1~]-\-q8 q(p — l) + pu 

u(x) < cd(x,dn)- » , v(x) < cd{ x , an) * — . 

Concerning system (M), we show that v always satisfies Harnack inequality. In the case 
SI = 5(0, 1)\{0}, we also study the behaviour near of the solutions of more general 
weighted systems, giving a priori estimates and removability results. Finally we prove the 
sharpness of the results. 
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1 Introduction 

In this article we study the nonnegative solutions of quasilinear systems in a domain Q of R , 
either with absorption terms, or mixed terms, that is, 

<"> (L1> 

where 5, ^ > and 1 < p, q < N. The operators are given in divergence form by 

A p u := div [A p (x, u, Vu)] , A q v := div [A q (x, v, Vv)] , 

where A p and A q are Caratheodory functions. In our main results, we suppose that A p is S-p-C 
(strongly-p-coercive) , that means 

A p (x,u,r]).7] > K ltP \r]\ p > K 2 , p \A p (x,u,r])\ p ' , V(x,u,rj) EfixK+x^. 

for some Kx )P , K2 tP > 0, and similarly for A q . The model type for A p is the p-Laplace operator 

u i — > A p u = div(\Vu\ p ~ 2 Vu). 

We prove a priori estimates of Keller- Osserman type for such operators, under a natural con- 
dition of " super linearity" : 

D = tf/i-(p-l)(g-l)>0, (1.2) 

and we deduce Liouville type results of nonexistence of entire solutions. We also study the 
behaviour near of nonnegative solutions of possibly weighted systems of the form 

(A \i ^p n = l x l a ' U<S > (M \ i A p u = \x\ a v & , 

1 W) \ A q V = \ X \ b , { W ' \ -A q V = \x\ b U» , 

in fi\ {0} , where a, b G R, with a > — p and b > —q. In particular we discuss about the Harnack 
inequality for u or v. 

Recall some classical results in the scalar case. For the model equation with an absorption 
term 

A p u = u Q , (1.3) 

in Q, with Q > p— 1, the first estimate was obtained by Keller |19| and Osserman [23] for p = 2, 
and extended to the case p ^ 2 in [29]: any nonnegative solution u£C 2 (fi) satisfies 

u(x) < Cd{x,dfl)- p/(Q - p+1 \ (1.4) 

where d(x,dQ) is the distance to the boundary, and C = C(N,p,Q). For the equation with a 
source term 

-A p u = u Q , 

up to now estimate (jl.4j) . valid for any Q > p — 1 in the radial case, has been obtained only 
for Q < Q* , where Q* = N ^~^ p +P is the Sobolev exponent, with difficult proofs, see [18], [9] 
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in the case p = 2 and [27] in the general case p > 1. For p = 2, the estimate, with a universal 
constant, is not true for Q = ^r|, and the problem is open between Q* and jf^g- 

Up to our knowledge all the known estimates for systems are related with systems for which 
some comparison properties hold, of competitive type, see [16], or of cooperative type, see [TT] : 
or with quasilinear operators in [IT], [32] , Problems (A) and (M) have been the object of very 
few works because such properties do not hold. The main ones concern systems (A w ) and (M w ) 
in the linear case p = q = 2, see [5] and [6]; the proofs rely on the inequalities satisfied by the 
mean values u and v on spheres of radius r, they cannot be extended to the quasilinear case. 
A radial study of system (^4) was introduced in [TB] , and recently in [7J . 



The problem with two source terms 

(5) 



-A p u = \x\ a v 5 , 

-A q V = 



was analyzed in [8j. The results are based on integral estimates, still valid under weaker 
assumptions: A p is called W-p-C (weakly-p-coercive) if 

A p (x,u,r]).r] > K p \A p (x,u,rj)\ p ' , V(x, u, rj) G x R+ x R N (1.5) 

for some K p > 0; similarly for A q . When 5, p < Q\, where Q\ = N ^Z p 1 ^ , punctual estimates 
were deduced for S-p-C, S-q-C operators and it was shown that u and v satisfy the Harnack 
inequality. 

In Section [2j we give our main tools for obtaining a priori estimates. First we show that 
the technique of integral estimates if fundamental, and can be used also for systems (^4) and 
(M). In Proposition 12 . 1 1 we consider both equations with absorption or source terms 

— ApU + / = 0, or — A p u = f, (1.6) 

in a domain Q., where / € Lj oc (Q), f > 0, and obtain local integral estimates of / with respect 
to u in a ball B(xq,p). When A p is S-p-C, they imply minorizations by the Wolf potential of 
/ in the ball 

i 

p-i 



w( p (B(x ,p))= f / ) y, (1-7) 

Jo \ JB(x a £) / * 



l(x ,t) 

extending the first results of [20], |21| . The second tool is the well known weak Harnack 
inequalities for solutions of (11.6P in case of S-p-C operators, and a more general version in case 
of equation with absorption, which appears to be very useful. The third one is a boostrap 
argument given in [5] which remains essential. 

In Section [3] we study both systems (A) and (M). When A p = A p and A q = A q , they 
admit particular radial solutions 

u*{x) = A* |x|" 7 ,v*(r) = B*\x\- ( , 

where 

= Kg -l)+qS = q{p -l)+pp 

7 D ' D ' 
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whenever 7 > ^— ^ and £ > for system (A), and 7 > and £ < f° r system (M). 

Our main result for the system with absorption term (A) extends precisely the Osserman- 
Keller estimate of the scalar case (11.31): 



Theorem 1.1 Assume that 

A p is S-p-C, A q is S-q-C, (1.8) 
and (TO)) ZioWs. Lei u € W^'c (^) H C (0) , v G (O) n C (U) be nonnegative solutions of 

-A p u + v s < 0, . _ 

-A^ + < 0, 

T/ien /or any xe!l 

u(x) < Cd(x,dn)-" / , v(x) < Cd(x,dn)~t, (1.9) 

with C = C(N,p,q,6,(j,,K l!P ,K 2> p,Ki )q ,K2,q). 

Our second result shows that the mixed system (M) also satisfies the Osserman-Keller 
estimate, without any restriction on 5 and [/,, and moreover the second function v always 
satisfies Harnack inequality: 

Theorem 1.2 Assume [IHj) . ITR) . Let u G Wj£ (O) D C (fi) , v G W^ 5 (O) n C (0) 6e nonneo- 
a£i?;e solutions of 

-A p u + v s <0, . _ 



T/ien siiZZ ZioZds /or any 

Moreover, if (u,v) is any nonnegative solution of system (M), then v satisfies Harnack 
inequality in £1, and there exists another C > as above, such that the punctual inequality holds 

(x) < Cv q - l (x)d(x, dn)- q . (1.10) 

Notice that the results are new even for p = q = 2. As a consequence we deduce Liouville 
properties: 

Corollary 1.3 Assume tl.ty ), [T78\) . Then there exist no entire nonnegative solutions of systems 
(A) or (M). 

Section [4] concerns the behaviour near of systems with possible weights (A w ) and (M w ), 
where 7, £ are replaced by 

_ (p + a)(g-l) + (q + b)5 _ (q + b)(p - 1) + (p + a)n 

in other terms o~£ a) 6 = (p — l)7o,6 + p + a, fi~/ a ,b = (l — l)£a,& + q + b. We set B r = B(0, r) and 
B' r = B r \ {0} for any r > 0. Our results extend and simplify the results of [5], [6] in a significant 
way: 
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Theorem 1.4 Assume fO )) .fO p . Zei u G W^(Si) n C , u € W/ G ' c 9 D C(Bj) &e 
nonnegative solutions of 

-A p u + \x\ a v s < 0, , . u 

A \ \h n ^ n 111 Bl. (1.11) 

Then there exists C = C(N,p, q, a, b, 6, fj,, Ki >p , K2 tP , Ki >q , if2,<j) > such that 

u(x) < C \x\~~ /a ' b , v(x) < C\x\^ a ' b inB[. (1.12) 

2 

Theorem 1.5 Assume fQj) .(TO j) . £e£ u G W/ 'f nC(5j), u G W^'J (-B() n C (B[) be 
nonnegative solutions of 

-A p u + \x\ a v 6 <0, . , 

in B[ . Then there exists C > as in theorem \l-4\ such that 

_ _c N ~i 

u(x)<C\x\ ~' a ' b , v(x) < C mm(\x\ u > b , \x\ i- 1 ), in B L . (1-14) 

2 

Moreover if (u,v) is any nonnegative solution of (M w ), then v satisfies Harnack inequality in 
B'i , and there exist another C > as above, such that 

2 

\x\ b+ V (x) < Cv q -\x), in B'i . (1.15) 

2 

Moreover we give removability results for the two systems (A w ) and (M w ), see Theorems 
14.11 14. 2^ whenever A p and A q satisfy monotonicity and homogeneity properties, extending to 
the quasilinear case [5j Corollary 1.2] and [BJ Theorem 1.1]. 

In Section [5] we show that our results on Harnack inequality are optimal, even in the radial 
case. And we prove the sharpness of the removability conditions. 

2 Main tools 

For any x G 1^ and r > 0, we set B(x, r) = G M / \y — x\ < r} and B r = B(0, r). For any 
function w G L 1 (0), and for any weight function ip G L°°(Cl) such that p > 0, ip ^ 0, we denote 
by 

(p w = -= — / wp 
the mean value of w with respect to p and by 

w = — — \ ir — h ir. 



o |0| jo ./ ! 



For any function g G L\ oc (VL), we shall say that a function u G W l( £ (U) satisfies 
— A p u ^ g in f2, (resp. ^, resp. =) 
if A p (x, u, Vu) G L^ oc (f2) and 

— / A p (x, u, Vn).V0 ^ / <7</), (resp. ^, resp. =) (2.1) 

for any nonnegative <fi G W 1 ' 00 ^) with compact support in fL 
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2.1 Integral estimates under weak conditions 

Next we prove integral inequalities for equations (|1.6p . either with source or with absorption 
terms. Such type of estimates for source terms were initiated in [8j Proposition 2.1]. They were 
used for specific / for proving Liouville type results in |23j and more recently in [10, Theorem 
2.1], and in p3] for the case of absorption terms. They are obtained by multiplication by u a , 
with a < for the source case, a > for the absorption case. 

Proposition 2.1 Let A p be W-p-C. Let f € L} oc (Vt), f > and let u £ (O) be any 

nonnegative solution of inequalaity 

- A p u ^ /, in O, (2.2) 



or 

- A p u + f S 0, in n. (2.3) 

Let £ £ T> (f2), with values in [0,1] , and ip = £ A , A > 0, and S% =supp\V^\. Then for any 
£ > p — 1, there exists \(p,£) such that for A > X(p,£), there exists C = C(N,p, K p ,£, A) > 
such that 

p-i 

J ftp S C|S 5 |max|V£| p \J) u^j . (2.4) 

Proof, (i) First assume that £ > p — 1 + a, with a £ (1 — p, 0) in case of equation (|2.2p . 
a G (0,1) (any a > if n S L^ C (Q)) in case of equation (|2,3p . We claim that there exists 
A(p, a, ^) such that for any A > \(p, a, £) 



p—l-\-g 

J fu a <p S C|S f |max|V£| p (j> u'tpj , (2.5) 



for some C = C(N,p,K p ,a,£, A). For proving (|2.5p . one can assume that u G L 1 (B(xq, p)). 
Let = £ A , where A > will be chosen after. Let 8 > 0, k > 1, and (r/ n ) be a sequence of 
mollifiers; we set u$ = u + 5, us t k = min(u, k) + 5 and approximate -u by us : k,n = ^<5,fc * "Hn, and 
we take </> = u^ k n p as a test function. Then in any case, from (II. 5h and Holder inequality, 



l a l / u 6kn ( P A p( x > u > S7u )- S7u 6An+ / f u 5,k,nf 

Jn ' ' Jn 
<A / < fcin £ A -V P (^,Vu)||V£| 
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Otherwise (Vus t k,n) tends to X{u<k}^ u m ^?oc(^)> anc ^ U P *° subsequence a.e. in Q, and 
A p (x,u,Vu) £ Lf oc (n). By letting n — > oo, we obtain 

M / u^ A A p (x,u,Vu).Vu + / /n^ A 
J{«<fc} ' 

i i 

with C = C(a, K p ,p, A); otherwise, for a < 1 (or it G L^ c (f2) and taking k > sup s u) 
/ u^ 1 £ x A p (x,u,Vu).Vu= / ^" 1 e A Ap(x,u,Vu).Vu+ / n^^ A A p (x, u, Vu).Vu 

Jfi ' J{u<k} ' J{u>k} ' 

^ / ^" 1 e A Ap(x,u,Vu).Vu + MA; a - 1 
where M = Jq £ A A p (x, u, Vu).Vu (or M = 0) is independent of k and (5. Then, for any 6 > 1, 
M /" 4-^ A A p (x, U ,V U ).V U + / /< fc £ A <C / u a + k p - l e- p \Vt\ p + Mk^ 



1 i 
6 ( r \W 



<5,fc 



a— 1 
p' 



Choosing 6> = ^/(a + p — 1) > 1, and A > A(p, a, £) = p6', we find 

M / ^- 1 £ A A p ( a ^,Vu).Vu + / /< fc £ A 
* J{u<k} Jn 

a-f p — 1 1 

< C U ' ^ jf |V^rj 6 + Mk^ 

a+p — 1 

< C|^max|V£| p (j u l 5 y\ +Mk^, 
with a new constant C = C(N,p, K, a, £). As k — > oo, we deduce 



a+p— 1 
1 



y^^-VA p (x,n,Vn).Vn + y^/n^<C|S f Kmax|Ve| p ^n^J . (2.6) 

Finally as 5 — > we get (|2.5p with a new constant C. Moreover we obtain an estimate of the 
gradient terms: 

( ct+p-l 
[ «V) 1 . (2.7) 
Jn J 



(ii) Next we only assume that t > p — 1, u G L 1 {B{xq, p)). Let as above, and fix 
some a = a(p,£) such that a € (1 — p, 0) and (1 — a){p — 1) < £ for ()2.2|) . a € (0,1) and 
a+p- 1 < £ for ([23]). In any case r = ^/(l - a)(p - 1) > 1, and l/Op' + l/pr = (p-l)/£. Let 
A > X(p, a(p, £), £) > pr'. We take if as a test function and from (|2.6p we deduce successively, 
with new constants C, 



ct — l 1 — g 



/ JV<a/ C A - X |A p (x, n,Vn) I I V^l <C / ^jA^z, «, Vu)| | V£| V 

i j_ 

<cn <-va p (o;,«,v«).v«j ^ n^" ^ e A - pr 'iver' 

< C \S^ + ^ m&x\VC\ p I J u\ 



4tt + — 

p'# pr 



p-i 



and dSH) follows as (5 0. ■ 

Corollary 2.2 Under the assumptions of Proposition \ 2.1l consider any ball B(xo,2p) C £1, 
and any e € (0, ^] . Let cp = £ A u>i£/i £ swc/i that 

£ = 1 inB(x ,p), 5 = 0«% +£) , |V£| < — . (2.8) 

£p 

Then for any £ > p — 1, i/iere exists A(p,^) such that for A > X(p,£), there exists C = 
C(N,p,K,£,X) > such that 

p-i 

f<C{e P yv^j>u^ 1 . (2.9) 

Remark 2.3 If = ujL 1 S'| where the are 2 by 2 disjoint, then (^\) can be replaced by 

k p- 1 
/ f<P S CV|5i|max|Ve| p | I u l \ . (2.10) 

2.2 Punctual estimates under strong conditions 

When Ap is S-p-C, the estimate (|2.7p of the gradient is the beginning of the proof of the 
well-known weak Harnack inequalities: 
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Theorem 2.4 (|25j, [28| ) Let A p be S-p-C, and u € (O) 6e nonnegative, such that 

—A p u ^ in f2; iaen /or any 6a// B(xq, 3p) C O, and any £ > p — 1, 



sup u < C * mM , (2.11) 

B(z ,p) V^B(xo,2p) / 



wni/i C = C(N,p, £, Ki >p , K2 iP ). Let w € (O) 6e nonnegative, such that —A p w > 0; ia 
/or any te// 5(xo,3p) C Jl, /or any £ € (0, iV(p - l)/(N-p)) 



CD 



v K \ <C inf u. (2.12) 

B(x ,2p) I B{x ,p) 

Next we give a more precise version of weak Harnack inequality (12. lip . Such a kind of 
inequality was first established in the parabolic case in |12j . 



Lemma 2.5 Let A p be S-p-C, and u E (f2) 6e nonnegative, such that —A p u ^ 0. Then /< 
any s > 0, taere exists a constant C = C(N,p, s, Ki }P , K2 }P ), such that for any ball B(x$,2p) C 
SI and any £ € (0, |] , 



or 



i 

\ s 



sup u < Ce~^~ <f> n s . (2.13) 

B(ko,p) \iB(x ,p(l+£)) / 

Proof. From a slight adaptation of the usual case where e = \, for any I > p — 1, there exists 
C = C(JV, ^) > such that for any e € (0, |) , 

sup u S Ce~ N [I u l J . (2.14) 

B(x ,p) \JB(io,(i(14«)) / 

Thus we can assume s < p — 1. We fix for example I = p, and define a sequence (p n ) by po = p, 
and p n = p(l + | + ... + (§) n ) for any n > 1, and we set M n = swp B r xo pn ^ u p . From (|2.14p we 
obtain, with new constants C = C(N,p), 

Mn ^ C(^±i - 1)-^ / uP< C(-)- (n+1)iVp / « p - 

P»» J B(x , Pn+ i) 2 J B{xo,pn+i) 

From the Young inequality, for any 5 £ (0, 1), and any r < 1, we obtain 



M n s cc-r {n+1)Np M^~i 



B(xo,p„+i) 



S SM n+1 + rd 1 - 1 / r '(C(-)- (n+1)7Vp )? 

2 B(xo,p n+ i) 

Defining « = r^-V'Cr and 6 = (f )- jv pA, we find 

M n ^5M n+1 + b n+l K\l u pr \ 

\J B(x ,p„ + i) I 



U 
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Taking 5 = ^ and iterating, we obtain 

M = sup vP S S n+l M n+1 + bnJ2(5by I <f 

B(x ,p) i=Q yJ_B(xo,Pn+i 



) 



1 



^ 5 n+1 M n+1 + 26k * u? r . 

\J B(x ,p n+ i) J 

Since B(xq, p n +i) C B(xq, p(l+s)), going to the limit as n — )• oo, and returning to u, we deduce 



sup u < {2bK) 1/p I I 

3(xn,p) \JB(xn,p[l- 



j_ 

rp 



B(x ,p) \J B(x ,p(l+e)) J 

and the conclusion follows by taking r = s/p. ■ 

It is interesting to make the link between Proposition ^.!) with the powerful estimates issued 
from the potential theory, involving Wolf potentials, proved in [20], [21] and [22]. Here we show 
that the lower estimates hold for any S-p-C operator. 



Corollary 2.6 Suppose that A p is S-p-C. If u is a nonnegative solution of 112. 2\) . then for any 
ball B(x , 2p) C f2, 

CWL(B(x ,p))+ inf u< liminfu(x), (2.15) 

' p B(x ,2p) x^-xo 



where W{ p is the Wolf potential off defined at fli. 7\ ), andC = C(N,p, K\ jP , i^2,p)- If u satisfies 
JO) , then 

CW( (B(xo,p)) + limsupu(x)< sup u. (2-16) 

x^tx B(x ,2p) 

Proof, (i) The function w = u— m 2p , where m p = mf B ( Xo p j u, is nonnegative in B(xo,2p), 
and satisfies the inequality —B p w > /, where 

w i — > B p w = div Ap(x, w + m2 P , Vw) 
is also a S-p-C operator. Then from Proposition 12.11 with £ as in (12 .8p . fixing ^ € (0, — - — — 



N-p 

and e = i, and applying Harnack inequality (I2.12p . there exists C = C(N,p,Ki :P ,K2 :P ) such 
that | | 

2C[p l - N [ f) <p- X f/ («-m 2p )M <p- 1 (m p -m 2p ). 

\ JB(xo, P ) J \Jb(xo,2 P ) J 

Setting = 2 1 - J ' / 9, as in [20], 



CW-j fJB(xo, p)) < y^(w , — m - ,) = limm , — inf u = liminf u — inf u. 

,P ' ' 3 B(x ,2p) a>-*Bo B(x ,2p) 

(ii) The function y = M 2p — it where M 2p = sup B ( a . 0j2p ) u satisfies the inequality —C p w > f in 
B(xq, 2p), where 

w i — C p it> := div [Ap(x, M 2p — w, Vw)] 
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is still S-p-C. Then 

W( (B(xq, p) S C{ sup u — limsup-u) 

B(x Q ,2p) X-Kco 

and ([2.160 follows. ■ 



Remark 2.7 The minorizations by Wolf potentials i2.15\) and \2. 1 6\) have been proved in 120^ 
and f2~m for S-p-C operators of type A p u := div [A p (x, Vu)] independent of u, satisfying more- 
over monotonicity and homogeneity properties, in particular A p (— u) = —A p u. The solutions 
are defined in the sense of potential theory, and may not belong to W,'^ (fi) , / can be a Radon 
measure; majorizations by Wolf potentials are also given, with weighted operators, see f21f and 
jl22f . In the same way Proposition \2.1\ can also be extended to weighted operators, see J3, Re- 
mark 2.4] and llJfl , or to the case of a Radon measure when A p is S-p-C by using the notion 
of local renormalized solution introduced in f3$[. 

2.3 A bootstrap result 

Finally we give a variant of a result of Lemma 2.2]: 



Lemma 2.8 Let d,h € R with d G (0, 1) and y, <3? be two positive functions on some interval 
(0, R] , and y is nondecreasing. Assume that there exist some K,M > and Eq G (0, J such 
that, for any e € (0,£o], 



yip) S Ke- h ${p)y d [p{\ + e)] and max $(r) <s M $(p), 

re[p,3f] 

Then there exists C = C(K,M,d,h,£o) > such that 



!j(i>)<('<lnn) — . , (0,^ 



v> e [ 0, 1 



(2.17) 



Proof. Let e m = e /2 m (m € N), and P m = (1 + ei)..(l + e m ). Then (P m ) has a finite limit 
P > 0, and more precisely P < e 2e ° < e. For any p E (0, ^] and any m > 1, 

y(A-i) < Ke7 n ^{pP m -i)y d (pP m ). 
By induction, for any m > 1, 

Hence from the assumption on <I>, 

and y dm (pP m ) < y dm (ep) < y dm (§), and ]imy dm (^) = 1, because d < 1. Hence ([2.170 follows 
withC = (^ /l ) 1 /(i-^) 2 V(i^) 2 M d /( 1 - rf ) 2 . ■ 
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3 Keller-Osserman estimates 



3.1 The scalar case 

First consider the solutions of inequality 

-A p u + cu Q <0, in ft, (3.1) 

with Q > p — 1 and c > 0. From the integral estimates of Proposition 12 . 1 1 we get easily Keller- 
Osserman estimates in the scalar case of the equation with absorption, without any hypothesis 
of monotonicity on the operator: 

Proposition 3.1 Let Q > p - 1, c > 0. If A p is S-p-C, and u G (ft) Pi C (ft) is a 

nonnegative solution of there exists with C = C(N,p, K\ tP , K2, p , Q) such that, for any 

x G ft, 

u{x) < Cc' 1 ^ Q+1 - p '>d(x,dn)- p ^ Q+1 -P\ (3.2) 

Proof. Let B(x ,po) C ft, and u G W x ' p (B(x ,p )) . From Corollary EjZ] with p < e = \, 
and £ = Q and a function (p satisfying (j2.8|) . we obtain for A = X(p, Q) 

\ — 

u Q <c- x Cp- p (J> u Q j Q , (3.3) 
where C = C(N,p, K\ tP , K2 iP , Q). Then with another C > as above, 

\ * i p 

u J <Cc Q+i-pp Q+i-p. 

Since A p is S-p-C, from the weak Harnack inequality ([2. lip , with another constant C as above, 

u(x )<C\(f u Q ] < c~ o+i-f p~ Q+i-p , 
\Jb(x ,p) J 

and (|3.2p follows by taking po = d(xo,dfl). ■ 

3.2 The system (A) 

Proof of Theorem ll.il Consider a ball B(xq, po) C ft, e G (0, |] , and a function ip satisfying 
(|2.8p with A large enough. 



(i) Case p > p — 1, 8 > q — 1. Here C denotes different constants which only depend on 
N,p,q,5, p, and K\ jP , i^2,p, -^2,q- We take e = | and apply Proposition 12.11 with p < ^ to 
the solution u with / = v s , and with £ = p > p — 1; since .Ap is W-p-C, from (|2.9p . we obtain 

p-i 

< ( I uA " , (3.4) 
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and similarly we apply it to the solution v with now / = and i = 5 > q — 1 : since A q is 
W-q-C, we obtain 



,5 X 



< Cp- q [<p v 6 ) . (3.5) 



We can assume that § u 1 * > 0. Indeed if § u* 1 = 0, then u = in B(xq,pq). Then Vu = 0, 
thus v s = and then the estimates are trivially verified. Replacing (|3.5p in (|3.4p we deduce 



« 5 < Cp- p - qE ^ (I v s 



(g-i)( P -i) 



and similarly for u, hence 



i , _ . i 



< Cp~t, [<p < Cp~ 7 . (3.6) 



Moreover, since *4. g is S-g-C, then from the usual weak Harnack inequality, since v 6 (J)), 
and </j(:e) = 1 in B(xq,p), with values in [0, 1] , 



sup v<CU v 5 < * V s ) < Cp-t. 

B(xo4) \Jb(x ,p) J \J<p J 

Similarly for u, since A p is S-p-C. 

(ii) Case p > p — 1, and (5 < q — 1. Here we still apply Proposition 12.11 with p < 
e G (0, 1/4] , and a function <p satisfying (|2.8p . Since /i > p — 1, we still obtain (|3.4p ; and for 
any /c > g — 1, and A large enough, 

/<%r , (3.7) 

and from Lemma [H 



v fe V < sup v<Ce'^(j) V s ) . 

lip J B(x ,p(l+e)) \J B{xoAWe)) J 

Then with new constants C, setting m = q + 5~ 2 Nq 2 (q — 1), and h = {p — l)p^ 1 m, 

/ \ ^ 

u p < Ce- m p- q [ I V s ) , (3.8) 

\J B(x ,p(l+2e)) J 

hence from (13.41) and (13.81). 



J v s < C 6 v 5 < Cp- p (in**) " < Ce~ h p~ P ^ P " f i / 

B(xo,p) Jf \Jf> J \J B{x Q ,p{l+2e)) 
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for any p < ^. Next we apply the boostrap Lemma [2.81 with R = po, y(p) = § B r XQ p \ v s , 



Pfi + q(p-l) 

<£(r) = r m and 2e. We deduce that 



v 1/6 
v 5 < Op'*, 



B(xo,p) 



for any p < ^e, and thus also 



sup v < C ( & v s \ < Op'*, sup «<CU uM < Cp~\ 



B(x 0! £) \JB(xo,p) 
In particular 

u(x ) < Cp-\ v(x ) < Cp^, (3.9) 
for any ball B(xq,Pq) C ft, and the estimates (|1.9|) follow by taking po = d(xo,d£l). ■ 

3.3 The system (M) 

Proof of Theorem 11.21 We consider a ball B{xq,pq) such that B(xQ,2po) C SI. From Propo- 
sition 12. H we have the same estimates: for any £ > p — l,k > q — 1, p < po, 

( 9-1 „ P-1 

<f v k ) " , <f v s <Cp- p (<f u e 

From Lemma 12.51 (even if p < p — 1), we have 

sup u" < C (f u*. 

B(x ,§) JB(x ,p) 

Taking k < N ^_^ , and using the weak Harnack inequality for v, we obtain 
sup /<Ci u* <C i u* < Cp- q ( & v k 1 

B(xo,f) JB(x ,p) J<p \J<p 

q-l 

<Cp- q \(f v k ) <Cp- q inf ufe" 1 ); 

\JB(x ,2p) J B(x ,p) 

hence (jl.lUp holds in B(xq, £). Moreover if v(xo) = 0, then u = in -B(xo, f)' then also v = 
in -B(xo, f )• Since ft is connected, it implies that v = 0, and then u = 0. If v ^ 0, then v stays 
positive in ft, and we can write 

-A q v = dv 9 ~ 1 , in ft, (3.10) 
with d(x) = u^/v^^ 1 ^ < Cp~ q in B(xq, ^); in particular 
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thus (jl.lOp holds and v satisfies Harnack inequality in 0, : there exists a constant C > such 
that 



Therefore 



sup v < C inf v. 

B(x ,p) B(x ,p) 



v 5 (x ) < sup v s < C inf v s < C <b v s < Cp- p <k u 



B(x ,p) B(x ,p) 



ip \Jtp 



p-1 



i „p-i ( g -i)(p-i) 

< Cp^ p sup u^ 1 < C p -' p p q ~ inf u i 

B(x ,2p) B{x 0) ip) 

KCp-^^v" 3 ^^); (3-12) 
and then (|3.9p follows again from f)3. 12|) and (|3.1ip . ■ 

Remark 3.2 Once we have proved h3.11\) we can obtain the estimate on u in another way: we 
have the relation in the ball 

A p u = v > cui- 1 , inB{xQ,po) 

— s 
with c = C\Pq 1 ; then from Osserman-Keller estimates of Proposition [3~l\ with Q = > p— 1, 

we deduce that 

u{x) < Cic-Wp' 75 ^ = C 3 Po 7 , in B(x , P j-). 

The Liouville results are a direct consequence of the estimates: 

Proof of Corollary 11.31 Let x € M, N be arbitrary. Applying the estimates in a ball B(x,R), 
we deduce that u(x) < CRT 1 ,v(x) < CR~^. Then we get u{x) = v(x) = by making R tend 
to oo. ■ 

Remark 3.3 In the scalar case of inequality h3. 1\) it was proved in \14f that the Liouville result 
is also valid for a W-p-C operator. In the case of systems (A) or (M), the question is open. 
Indeed the method is based on the multiplication of the inequality by u a with a large enough, 
and cannot be extended to the system. 

4 Behaviour near an isolated point 

4.1 The system (A w ) 

Proof of Theorem II. 4L It is a variant of Theorem ll.il we consider Vt = B[ and xo £ B' x , and 

2 

take po = -^p. Here we apply Proposition 12.11 in the ball B(xo,p) with p < ^ and e € (0, |] . 
The estimates (|3,4|) and (|3.7|) are replaced by 

\x\ a v s < C(ep)- p ( I u l j 1 , i \x\ b u^ <C(ep)- q ( (f v k ) * , (4.1) 
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for any £ > p — l,k > q — 1; and 2po < < 6po in B(xq, 2po), then in any of the cases a < 
or a > 0, with a new constant C, 



/ < C£- p /T (p+a) ( / i/V , / < Ce-V _(9+b) ( / vA " . (4.2) 

Then all the proof is the same up to the change from p, q into p + a and q + b. We deduce the 
same estimates with 7, £ replaced by 7 a ,&,£a,6 : 

-u(xo) < C|x |" 7a ' 6 , u(x ) < C\x \~ $a ' b , (4.3) 
where C depends on N,p, q, a, b, 5, p, and -Ki, p , AT2 iP , Ki,q,K2,q- ■ 

4.2 The system (M ro ) 

Proof of theorem II. 5L In the same way we obtain estimate (14.3p . then we only need to prove 

N-q 

the estimate with respect to \x\ 9-1 . We can apply to the function v the results of [2], recalled 

in P Propositions 2.2 and 2.3]: G L 1 fBi) , and for any k G (0, , and p > 
small enough, 

(/ u fe ] <Cp~^. (4.4) 
\Jb(o,p) J 

Moreover, arguing as in the proof of ([l.lOp . we obtain the punctual inequality 

u^(x ) < C \x r {q+b) v^ixo), in B[, (4.5) 



a(sco) = Fo -777 — r < C Fo y . 



which implies that 
Then v satisfies the Harnack inequality in B\ , hence from (|4.4p 



v{x ) < If tr < C \x Q \ 1- 1 



and (fLTil follows. 



4.3 Removability results 

Here we suppose that 

{A p u := di-u [A p (x, Vu)] , A p is S-p-C, 

(A p (z,0-A p (x,0).(£-0 >0, for^C, 

Ap(x, AC) = |A| P ~ 2 AA p (x, 0, for A / 0, 

and similarly for We give sufficient conditions ensuring that at least one of the functions 
u, v or both are bounded. We have the two following results, relative to systems (A w ) and 

(M w y. 
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Theorem 4.1 Assume IFE . (C p ),(C q ). Let u G (B[) , v G (B[) be nonnegative 

solutions of 

-ApU+ \x\ a v s < 0, . , 

-A q v + \x\ b u» <0, m v 

(i) If ^ a ,b < then u is bounded near 0; if Ca,b < then v is bounded. 

(ii) If moreover (u,v) is a solution of (A w ) and u is bounded near and 5 > fr+ff - 1 ) (or 
$ = (g+^li. 1 ) = A,pJ i/ien v is a/so bounded. In the same way if v is bounded and 
A* > ^ (or ii= ^ q+ N^£ p ^ if A q = A q ) then u is also bounded. 

Theorem 4.2 Assume (Tty . (C p ), (Q). Let u G Wj£ (B[) n C (B£) , u G Ty/ ' c 9 n C (B[) 
be nonnegative solutions of 

-AnU + \x\ a v s < 0, . „, 

-A q V > \x\ b U^ , 1 

Ifla,b < frf > or */7a,fe > frr anrf M > ^7^—^ » ^en u is bounded. 

The proofs require some lemmas, adapted to subsolutions of equation A p u = 0. 

Lemma 4.3 Assume (C p ). Let u G W l( £ (B[) n C(B[) 6e nonnegative, such that —A p u < in 
B[. Then either there exists C > and r G (0, |) suc/i i/iai 

p-JV 

sup it > Cpp- 1 , for any p G (0, r) , (4-6) 

\x\=p 

or u is bounded near 0. 

Proof. From our assumptions on A p , there exists at least a solution E of the Dirichlet problem 

— ApE = So, in Bi, 

where 5q is the Dirac mass at 0, in the renormalized sense, see |X3|, Theorem 3.1]. In particular 
it satisfies the equation in T>'(Bi), and it is a smooth solution of equation A P E = in B[. From 

N—p _ N—p 

[25] . [26] , there exists 61,62 > such that C\ \x\ p- 1 < E(x) ^ C2 |x[ p- 1 near 0. Assume 
that (I4.6P does not hold. Then there exists r n < min(l/n, r n _i) such that 

I p~ N \ 
sup u < -r n p_1 < ——E(r n ). 
\x\=r„ n nUi 

Next we use the comparison theorem in the annulus C n = jx G : r n < \x\ < ^} for functions 
in Wl£ (C) n C((Q, and we find that 

itfsc) < —prB(x) + max u, in C n . 

«6i |ar|=-| 

Going to the limit as n — > 00, we deduce that u is bounded. ■ 

Our next lemma complements the results of Proposition 2.2]: 
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Lemma 4.4 Assume that A p is W-p-C. Let f G Lj oc (B[), f^O.Letue Wj£(B[) be nonneg- 
ative, such that 

-A p u + fSO inB[. 

N-p 

If \x\p~ 1 u is bounded near 0, then f G L\JB X ). 

Proof. Let < p < \. Here we apply Proposition 12. II with ip = £ A given by 

1 3 C 

f = 1 for p < \x\ < — , f = for \x\ < —or Ixl > — , |Vf I < — . 

2 1 1 - 2 1 1 - 4 p 

From Remark 12.31 we find with for example i = p, 

p — 1 p — 1 

/ f^Cp N -4l J] 1 +c(l J] 1 . (4.7) 
Hence from our assumption on u, the integral is bounded, then / G L x (i3i). ■ 



2 

■V-p 



Proof of Theorem 14.11 (i) Suppose that 7^ < p ~f ■ Then u(xq) < C \xq\ p-* . Let us show 
that u is bounded. If ~f a ,b < it is a direct consequence of Lemma l4~3l Then we can assume 
la, b = ^rf • 11 u 1S n °t bounded, then f)4.6j) holds for some C > 0. Let us set / = |x| a u 5 . From 
(|4.2p with e = 4 then for any tq < \ and any xq such that \xq\ = vq, and Lemma 12.51 taking 



a = ro 



u"(*o) <cl u»< Cr {(!+b) - N ^ ( [ 

Jb(x , P ) \Jb(xo,2 P ) 



2-1 

5 



<Cr {q+b) - {N+a) ^[ I ^, / 



g-l 
S 



then 



Cr ^ b = Cr {q - 1] ^ b - q ~ b < sup u» < 0^+^+^ 

M=ro 

Cr J ' 1 J = Crl = C < I _ /; 



2-1 
5 



then for any n G N, 

C < 

By summation it contradicts Lemma 14.41 Similarly for v. 



(ii) Suppose that (u,v) is a solution of (A w ) and u is bounded and 5 > 1 • Here t> 

satisfies equation .4 g ?; = 5 with g = \x\ b u M ^ C , thus g G L Ar / <?+£ (r2) for some e > 0, then 
from [25J, [26J, if v is not bounded near 0, then there exist C\, C2 > such that 

N-g N-q 

Ci \x\ 1- 1 < v < Co \x\ 9-1 
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near 0. If S > {p+ f q - 1] then 

sr X N-q 

A p u = \x\ a v 5 > Ci\x\ a ~ b — = d\x\- p - £ , 
for some e > 0, then from JO}, 

p-i 

p- p - £ ^ C j> \x\- p - £ < Cp~ p (J> u^j 1 ^Cp- p , 

which is a contradiction. If 5 = ^ p+ ^j_~ , then 

C 2 \x\~ p > A p u = \x\ a v 5 > C!\x\- p . 

Otherwise u is bounded by some M in a ball B' r . Then the function w = M — u is nonnegative 
and bounded and satisfies 

-A p w>Ci\x\~ p in B' r . 

But for A p = A p , there is no bounded solution of this inequality, from [8j Proposition 2.7], we 
reach a contradiction. ■ 

Remark 4.5 The results obviously apply to the scalar case, finding again and improving a 
result of Jgl]/. 

Proof of Theorem 14.21 (i) Assume 7^ < ^E±- The proof of part (i) of Theorem 14.11 is still 
valid and shows that u is bounded. 

(ii) Assume 7^ > ^5f and p > ^ N+ N- P p ^ • Then £ a) 5 > ^Ei> thus the estimate (|1.14p for v 

N-q 

gives v(xq) < C \xq\ 9-1 , then 

u"(x ) < C \x \- {q+b) v^ixo) < C|x r (7V+fc) . 

N-p 

Then pp- 1 sup^^u tends to 0, hence u is bounded from Lemma 14.31 

Remark 4.6 Let us give an alternative proof of (i): the punctual inequality |^.5p implies that 
near 0, 

A p u > \x\ a v s > C\x\ a+ ^ q+b ^^u^^ q -^- 
then we are reduced to a simple scalar inequality: 

-A p u+ \x\ m u Q < 0, (4.8) 

with Q = ^tj- > p — 1 and m = a + -^pj^ > —p- And 7^5 = < ^E±l applying Theorem 

\4-l\ to the scalar inequality \4.8\) , we find again that u is bounded. 
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5 Sharpness of the results 



In this last section we show the optimality of our results by constructing some radial solutions 
of systems (A w ) or (M w ) in case A p = A p ,A q = A q . They are based on the transformation 
introduced in [3], valid for systems with any sign: 

-A p u = — div(|Vu| p_2 Vit) = E\ |x| a v s , 
—A q v = - div(| Vv\ q ~ 2 Vu) = e 2 \x\ b u'\ 

with Ei = — 1 = £2 for the system with absorption, and £\ = — 1, £2 = 1 for the mixed system: 
setting 



X(t) = , Y(t) = , Z(t) 

u v 



where t = lnr, and we obtain the system 



(EM 



x t = x 

Y t = Y 



X 
Y - 



p-1 " r p-1 

Af-g , _W_ 
g-1 ~T g-1 



Z t = Z [N + a - 5y - Z] , 



And u, v are recovered from X, Y, Z, W by the relations 



u- 



••'(IXI*- 1 Z)' 1 / D (\Y\ q - 1 W) 



{p-1) /D 



(5.1) 



5.1 About Harnack inequality 

Here we show that Harnack inequality can be false in case of system (A w ) and also for the 
function u of system (M w ), even in the radial case; indeed we construct nonnegative radial 
solutions of system (A w ) in a ball such that such that n(0) = < v(0), or by symmetry 
u(0) > = v(Q) and solutions of system (M w ) such that u(Q) = < v(0). Such solutions 
were constructed in [15] by using Schauder theorem, and in [7j in the case of system (Aw) for 
p = q = 2 by using system (£). Here we show that the construction of [7] extends to the general 
case. We consider the radial regular solutions, which are C 2 if a, b > 0, and C 1 if a, b > —1. 

Proposition 5.1 A p = A p ,A q = A q . For any vq > 0, there exists a regular radial solution of 
(A w ) and (M w ) such that u(Q) = < v(0) = vo. 

Proof. The regular solutions (u, v) with nonnegative initial data (uq, Vq) 7^ (0, 0) are increasing 
for system (A w ), hence X, Y < < Z, W and u is increasing and v is decreasing for system 
(M w ), hence X < < Y and Z, W > 0. As shown in [3J, the solutions (u,v) with it(0) = uo > 
0, u(0) = f > correspond to the trajectories of system (S) converging to the fixed point iVo = 
(0, 0, N + a, N + 6) as i — >• — 00, and local existence and uniqueness holds as in [JJ Proposition 
4.4]. As in [7] the solutions such that uq = < vq correspond to a trajectory converging to the 
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point S = (X, 0, Z, W) = ( — pizi, 0, N + a, N + b + p^j . The linearization at S gives the 
eigenvalues 

X 1 = X<0, A 2 = —— (q + b + p 1 -^-) > 0, A 3 = -^<0, A 4 = -TV < 0. 
q — 1 p — 1 

Then the unstable manifold V u has dimension 1 and V u fl{Y = 0} = 0, thus there exists a unique 
trajectory such that Y < (resp. 1" > 0) and Z, W > 0. There holds limt^-co e -A2 *Y = c > 0, 
limX = X, limZ = Z, limTV = W, then from (|5.ip v has a positive limit vq, an d u tends to 
0. By scaling we obtain the existence and uniqueness of solutions for any vq > 0. ■ 



5.2 About removability 

Here also we show that the results of Theorems 14.11 and 14.21 are optimal, by constructing 
singular solutions when the assumptions are not satisfied. We begin by system (A w ), extending 

N-p 
p-1 



Proposition 3.2]. Obviously it admits a particular singular solution when 7^ > ^— P and 



£a,& > — Moreover we find other types of singular solutions: 
Proposition 5.2 Consider system (A w ) with A p = A p ,A q = A q . 



0. 



(i) If p < — -, there exist solutions such that lim p _j.o P p_1 u = a > 0, lim p ^o v = f3 > 

(ii) If 5 < ^-^^j ^ and < ^ N+ N^ p p ^ , there exist solutions such that lim p _ i .oP p_1 u = 

N-q 

a > 0, lim^o p q ~ 1 v = (3 > 0. 

(in) If "f a ,b > Trrfi an d either p > ^"ffii^ ^ or p < 1 \ there exist solutions such 

N—p 1 / N—p _/ 1 t\\ 

i&ai lim p ^ p p - 1 « = a>0, lim p ^ P 9 " 1 1 ;J f = P{a) > 0. 

TTie results extend by symmetry, after exchanging u, v,a,j a b and v,u,b,£ a ,b- 

Proof. As in [5], [7] we prove the existence of trajectories of system (£) and return to u,v by 
using (|5.ip . 

(i) Such solutions correspond to trajectories converging to the fixed point G = (|rf , 0, 0, iV + 
b — ^Ef-p) °f (S). The linearization at Go gives the eigenvalues 

N-p , 1 , A^-p . A , iV-p Ar , 

Ai = f > 0, A 2 = -(g + 6 rM), A3 = TV + a > 0, A 4 = £/i — N — b. 

p — 1 g — 1 p — 1 p — 1 

If A» < {q+ N- P ~ 1] , then A 2, A 4 < 0. Then V„ has dimension 3, and V u n{Y = 0} and V u n{Z = 0} 
have dimension 2. This implies that V u must contain trajectories such that Y, Z < < X, W. 

(ii) Such solutions correspond to the fixed point Aq = f^5f , ^q~ii 0' • All the eigenvalues are 
positive: 

N-p N _ q N -q N-p 

Ai = -, A 2 = -, \ 3 = N + a- 6 -, A 4 = N + b - p -. 

p — 1 q — 1 q — 1 p — 1 
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The unstable manifold V u has dimension 4, then there exists an infinity of trajectories converging 
to A with X;Y,Z,W < 0. 

(hi) Such solutions correspond to the fixed point Pq = (^Ei, Y*, 0, W*^ = (jj^Ei,Y*,0,W*j , 
with Y* = - (g + b)), W* = N + b- fff/x. the eigenvalues are 

N - v D N - v 

F >0, A 2 = n, A 3 = -(7 f)>0, A 4 = -W*. 



p — 1 q — 1 p — 1 

-1 



If > "ffi^l 1 , then A 2 ,A4 > and thus V u has dimension 4, then there exist trajectories, 



with X,Y,Z,W < 0, converging to Po. If A* < ^tv-^T^ > then A 2 ,A4 < 0, V M has dimension 
2, and V u n {Z = 0} has dimension 1, thus there also exist trajectories with X, Z, W < < Y 
converging to Po- ■ 

In the same way, system (M w ) has a particular singular solution when j aj b > and 
ia,b < ^rf) an d we find other singular solutions: 

Proposition 5.3 Consider system (M w ) with A p = A p ,A q = A q . 

jy N — q 

(i) If^ia,b > > and ^ a) b > -q~ri> there exist solutions such that lim p ^oP q ~ 1 v = (3 > and 
limp^o p^ { ^ S ~ (q+a)) u = /3(a) > 0. 

(ii) If 5 < — j^z^ ^ and p < ^—7^^—^, there exist solutions such that limp_s.oP p_1 u 

X ; 

a > 0, lmip^o P 9-1 v = (3 > 0. 

Proof, (i) These solutions correspond to the fixed point Qo deduced from Po by symmetry, and 
our assumptions imply 5 > ^ N+ ^ 9 q ^ , hence there exist trajectories, such that X, Y, Z < < W 
converging to Qo- 

(ii) The conclusion follows as in Proposition 15.21 (ii). ■ 

We refer to [5] and [6] for a description of all the (various) possible behaviours of the solutions 
in the case p = q = 2. 
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